Quantum aging in mean-field models 
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We study the real-time dynamics of quantum models with long-range interactions coupled to a 
heat-bath within the closed-time path-integral formalism. We show that quantum fluctuations de- 
press the transition temperature. In the subcritical region there are two asymptotic time-regimes 
with (i) stationary, and (ii) slow aging dynamics. We extend the quantum fluctuation-dissipation 
theorem to the nonequilibrium case in a consistent way with the notion of an effective temperature 
that drives the system in the aging regime. The classical results are recovered for ft ^ 0. 
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The dynamics of nonequilibrium systems is being in- 
tensively studied now. Notably, glassy systems below 
their critical temperature have a very slow evolution with 
nonstationary dynamics [Q . Several theoretical ideas |^] 
are used to describe it, namely, scaling arguments, phase 
space models, analytical solutions to mean-field models 
and numerical simulations. The analysis of simple mean- 
field models (with long-range interactions) has provided 
a general scenario that is now being verified numer- 
ically for more realistic models All these studies 
concern classical systems. 

Recent experiments have motivated a renewed inter- 
est on the effect of quantum fluctuations (QF) on glassy 
systems. Up to present, theoretical studies have focused 
on how QF affect their equilibrium properties [p|~pd| . 

Since glasses below Tg are not expected to reach equi- 
librium in experimentally accessible times, it is important 
to device a method to understand the influence of QF on 
the truUy nonequilibrium real-time dynamics of this type 
of systems. Intuitively, one expects QF to only affect the 
short-time dynamics; however, they are also expected to 
act as thermal fluctuations. It is then not clear whether 
QF would destroy glassiness or modify it drastically. 

Our aims are (i) to present a formalism suited to study 
the real-time dynamics of a nonlinear, possibly disor- 
dered, model in contact with a bath; (ii) to propose a 
framework to study its dynamics that could be also ap- 
plicable to more realistic, finite dimensional, models; (iii) 
to show that below a critical line QF do not destroy the 
nonequilibrium effects of the glassy phase and that QF 
add up to an effective temperature Teff; (iv) to prove that 
Teff is non-zero even at zero bath temperature, that it 
drives the dynamics at late epochs and that it makes the 
dynamics appear classical in that time-regime. A longer 
account of our results will appear elsewhere. 

The real-time dynamics of a quantum system is de- 
scribed with a closed-time path-integral generating func- 
tional (CTP-GF) iQ. We choose a set of noninteracting 
harmonic oscillators with an adequate distribution of fre- 
quencies as a bath, and a linear interaction between bath 



and system The bath variables are next integrated 
out and the effect of the bath manifiests in the effective 
action through two non-local kernels associated to dissi- 
pation (?/) and noise (j/). If the model is disordered, one 
needs to compute averaged expectation values. However, 
the CTP-GF without sources is independent of the real- 
ization of disorder and one can hence avoid, in a quantum 
dynamical calculation, the introduction of replicas. In 
addition, when h ^ 0, the CTP-GF yields the classical 
Martin-Siggia-Rose one. 

For the sake of concreteness, we study the p spin-glass 



1 



N 



2mN 



En? 



E 

ii,...,ip 



(1) 



with Hi the canonical momenta, [11^,0^] = —ihSij. 
The multispin interactions Ji-^...ip are taken from a 
Gaussian distribution with zero mean and variance 
J2p!/(27VP-1) and EZiMit))]j = N, Vi. Square 
brackets denote an average over disorder and (•) the 
average over temporal histories. The quantum mean- 
field equations follow from a saddle-point approxi- 
mation and involve the symmetrized auto-correlation 
NC{t,tw) = [{(f){t)4>[tw) + 4>{tw)4>{i))]j and the response 
to an infinitesimal perturbation h applied at time t^, 
NR{t,t^) = S[((j){t))/dh{t^)]j\h=o. We define the ver- 
tex and self-energy: 

Dit, t^) = b{t, t^) + 2nv{t - t^) 
E(t, i™)-S(t, i^) - 477(t - t„) 

respectively. S and D are functions of C and R that fol- 
low either from the average over disorder or from approx- 
imations of the nonlinear interactions in the model (e.g. 
mode-couling approximations). For the p-spin model and 
t > t-u] they read 
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D{t: t^) + —^{t, t^) = ^— \ C{t, t^) + —R{t, t^) 
The terms associated to the bath are 
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1 



vit — tiu) = / dujl{u}) coth(/3?iijj/2) cos([x)(t — t^)) , 
Jo 



rj{t — tw) — 9{t ~ tw) f dijjl{ijj) sa\{ijj{t — t^)) . 
Jo 

We chose an Ohmic distribution of osciUator frequencies 
/(w) = {Mjo/tt) u) exp(— |w|/A) with A a cut-off and 
M70 a constant that plays the role of a friction coefficient. 
Their particular form depends on the choice of bath. The 
equations, for a random initial condition, are 

R{t,tw) = Go(t,tw) + / / dudvGo{t,u)T,{u,v)R{v,tjjj) 
Jo Jo 

C{t,t.u,)= [ du [ dvR{t,u)D{u,v)R{tni,v) (2) 
Jo Jo 

with G'^^{t,u) = 6{t - u) (md^ + ^(t)) the prop- 
agator and fi{t) either a Lagrange multiplier en- 
forcing a spherical constraint or the strength of a 
quadratic term. In the first case it is determined 
by the 'gap' equation: fi{t) — duYi{t,u)C{t,u) + 
D{t,u)R{t,u) — m 9j2C(t, tto)|j ,^t-- Causality implies 
R{t,t^) = ^{t,t^) = if L > t. The inertial 
term imposes continuity on the equal-times correlator 
Ihnt^^i- dtC{t,tw) = lhnf^^f+ dtC{t,tw) = 0, and 
R{t,t) = 0, \imt^^-,tdtR{t,tw) = 1/m. In what follows 
we set t > t^. The coupling to the bath implies dissipa- 
tion; if M70 ^ the energy density £ of the system de- 
creases. One can then envisage to switch-off the coupling 
(set M-fo = 0) when £ reaches a desired value and follow 
the subsequent evolution at constant £. This would be 
useful to further understand the energy landscape. Here, 
we keep Mjo 7^ for all times, reparametrize time ac- 
cording to t ^ M^ot and transform ?i in a free parameter 
n Mjoh. Consistently, C ^ C, R~^ R/{M-/o), m 
{M^oYm, J J, P P,A ^ A/(Af7o) and the units 
are [C] = [R] = [h] = 1, [m] = [/?] ^ [1/J] ^ [1/A] = [t]. 

We focus on p > 2, since p ~ 2 needs a special treat- 
ment Since Eqs. (^) are causal, we contruct its solu- 
tion numerically, step by step in t, spanning < t-^j < t. 
The numerical and analytical studies yield: 

(i) Quantum mode- coupling equations. For h and T 
above a critical line, there is a finite equilibration time t^q 
after which equilibrium dynamics sets in. The solution 
satifies invariance under time-translations (TTI) and the 
quantum fluctuation-dissipation theorem (QFDT). It is a 
'paramagnetic/liquid' phase. A TTI-QFDT ansatz yields 
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Riuj) = 

C{uj) = D{ijj)\R{oj)\'^ . 



(3) 

(4) 



with the couples R and C, and S and D, verifying QFDT: 
i?M = |^M/ge— tanh(^)c(.) (5) 



(r = t — tw)- Equations (||)-(0) are the quantum ver- 
sion of the mode-coupling equations used to describe 
supercooled-liquids [|^. Away from the critical line, C 
and R decay to zero very fast, with osicUations. Ap- 
proaching the critical line Tci{hd), the decay slows down 
and, if Td{h.d) 7^ 0, a plateau develops in C. At the 
critical line, the lenght of the plateau tends to infinity. 
The quantum critical point {Td — 0,hd ^ 0) is discussed 
below. 

(ii) Dynamics in the glassy phase. For h and T below 
the critical line, Teq — > 00 (as a function of N): times 
are always finite with respect to Teq. The system does 
not reach equilibrium. This is a 'spin-glass/glass' phase. 
There are two time regimes with different behaviours ac- 
cording to the relative value of i ~ and a characteristic 
(model-dependent) time T{t^): 

li t — tyj < T{tyj) {C{t,tyj) > q) thc dyuamlcs is sta- 
tionary; TTI and QFDT hold. In other words. 



g -f- CsT(t - tiL.) lim C{t,tm) 
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and RsT{t — t^) = limj^^oo Rit, tw)- The correlation ap- 
proaches a plateau q since limt_(^^oo CsT(i — tw) = 0. 
The response approaches zero, limt_t^_>oo ^st(^ — t^) — 
0. The equations for Cst and i?sT are identical to Eqs.(^)- 
(^) apart from contributions to jioo coming from the ag- 
ing regime. Cst and i?sT are linked through Eq.(||). 

lit—tw > T{tw) (C{t,tyj) < q) the dynamics is nonsta- 
tionary; TTI and FDT do not hold and there is quantum 
aging. The correlation decays from 5 to and we call 
it Ct^G{t,tyj). The decay of C becomes monotonous in 
the aging regime; the properties of monotonous two-time 
correlation functions derived in Ref. W imply: 
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with h{t) a growing function of time. For the p-spin 
model, is the identity. 

The system has a weak long-term memory, meaning 
that the response tends to zero but its integral over a 
growing time-interval is finite. The QFDT is more in- 
volved than the classical one since it is non-local (Eq.(||)). 
In the nonequilibrium case, we generalize it to 



R{t,ty,) = ^ 6{t-t^) 



du 
2^ 



exp(— ia;(i — tin)) 



C{t,Lu) 



xtanh('^(^4^^) Cit,u) (8) 



2Re / dsex.p{iuj(t - s))C{t,s) (9) 
'0 



with X{t, tin) a function of t and If the evolution is 
TTI and X{t,t^,) = 1 for aU times, Eq.(|) reduces to 
Eq.(||). Interestingly enough, Teff = T/X(t,tin) acts as 
an effective temperature in the system . For a model 
with two two-time sectors we propose 
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with a non-trivial function of h and T. When t 
and are widely separated, the integration over lu in 
Eq.(||) is dominated by w 0. Therefore, the factor 
ta,nh{Xj,c,{t,tw)l3hu}/2) can be substituted by X^QfihLo/2 
{even at T = if Xj^c{h,T) = x{h)T when T ~ 0). 
Hence, 



RAcit, tw) ^ (^{t — tw)Xf^c,f3dt^CAGitj tw) 



(10) 



and one recovers, in the aging regime, the classical mod- 
ified FDT Numerically (and experimentally), 
it is simpler to check it using the integrated response 
xit^tw) = J* dt"R{t,t") over a large time window 
that reads 

X[t,tw) ~ j J^-'dT'R,,{T') + ^iq^C,cXt,tw)) 

the first line holds for C{t,tw) > q while the second 
holds for C{t,tw) < q- Remarkably, the correlation and 
the violation of QFDT are given by similar expressions 
to the classical - though the values of q and Xag de- 
pend on fi. In this sense, we say that T(tw) acts a 
time-dependent 'decoherence time', beyond which the 
nonequilibrium regime is 'classical'. 




FIG. 1. The correlation function C{t + tw,tw) vs r for 
the p = 3 SG model, A = 5, J = l, m = l, r = 
and h — 0.1. The waiting times are, from bottom to top, 
— 2.5, 5, 10, 20, 40. q ~ 0.97. In the inset, the same curves 
for tw = 40 and, from top to bottom, h = 0.1, 0.5, 1, 2. 



An ansatz like Eq.(^ for Cag combined with Eq.(|T 
for the FDT violation solves the equations in the agin: 
regime. It yields 



Xao/T : 



{Rsriuj ^ 0))^ JMp ' 
i?sT(w = 0)(p-2)/g 
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(11) 
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that become J'^p{p — l)/2qP 
ip - 2)(1 - q)/q when ?i ^ §. 

QF depress the critical temperature. The transition 
line Td{hd) ends at a quantum critical point {T^ = 0,h(i ^ 
0). Equations (pT|)-(p^) indicate that when the transi- 
tion occurs at Td{'hd) 7^ it is first-order with Xj^q — > 1, 
Q ^ Qd and a finite linear stationary susceptibility 
i?sT(t^ = 0) < +00 (as in the classical limit). On this 
line, qd oc T^^^ and qd ^ for Td 0. At the quan- 
tum critical point, Eqs.(pT|)-(p^ suggest a second-order 
transition; if qd tends to zero as qd ~ {hd - ;i)"P/^ then 
Xac/T - {hd - ?i)-", and Rsrito = 0) ~ [hd - ;i)"(i-P/2) 
diverges when p > 2 (a is positive) . 

Let us now discuss the numerical results. In all figures 
p = 3, A = 5, J=l and ni — 1. We use T = and 
?i = 0.1 to illustrate the dynamics in the glassy phase. 
We also discuss the dependence upon T and h. 




FIG. 2. The response function for the same model as above. 
The waiting-times increase from top to bottom. In the in- 
set, check of FDT in the stationary regime. The full line is 
R{t,tw) for t = 40 fixed and tw € [0,40]. The dots are ob- 
tained from Eq.(^) with Xst = 1, using the numerical data 
for CsT(t - tw) = C{t,tw) - q (q ^ 0.97, see Fig.§). In both 
cases the response is plotted against t — tw- 
in Figs. 0and| we show the correlation C(t + tw,tw) 
(log-log plot) and the response R{T + tw,tw) (linear plot) 
vs the subsequent time r. These plots demonstrate 
the existence of the stationary and aging regimes. For 
t - tw < T{tw) (e.g. r(40) - 5) TTI and FDT are 
established while beyond T{tw) they break down. For 
h = 0.1 the plateau in C is at g ~ 0.97. C oscillates 
around q but is monotonous when it goes below it. In 
the inset we present the dependence of g on ?i for T = 0. 
QF generate a g < 1 such that the larger h the smaller 
q. The addition of thermal fluctuations has a similar ef- 
fect, the larger T, the smaller q. In order to check FDT 
in the stationary regime, in the inset of Fig. |^ we com- 
pare R{t,tw) from the numerical algorithm for t — AQ 
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fixed and e [0, 40] (full line) with R{t, t^) from Eq.(|) 
with X = I using CstH - tyj) = C{t,tw) - q, q ^ 0.97 
obtained from the algorithm (dots). The accord is very 
good at — tiu < T{tw) ~ 5. 

In Fig.^ we plot the integrated response x vs. C in a 
parametric plot. When C < q ^ 0.97 the x vs C curve 
approaches a straight line of finite slope I/Teff = X^q/T 
that we estimate from Eqs.(^ and ( p^ as Xj^^/T ~ 
0.60. The dots are a guide to the eye representing it. 
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FIG. 3. The integrated response xi"^ ^'tiii,tw) vs. the sym- 
metrized auto-correlation function C(T+tw,tw) in a paramet- 
ric plot. The curves correspond to = 10, 20, 30, 40. T = 
and h = 0.1. The dots are a guide to the eye and represent the 
analytic result Xf,c,/T ^ 0.60 for the slope of the asymptotic 
curve in the aging regime C < q ^ 0.97. 

To summarize, this formalism provides a clear frame- 
work to study the nonequilibrium regime and/or the 
eventual approach to equilibrium of quantum systems. 
It circumvents the difficulties of performing numerical 
simulations in real-time and predicts the existence and 
the structure of a rich nonequilibrium regime for glassy 
systems even in the presence of QF. Based on the suc- 
cess of mean-field like classical glassy models to describe 
some aspects of the dynamics of realistic glassy systems 
(see and references therein) one can expect that 

quantum mean-field models of the type considered here 
do also capture important features of the real-time dy- 
namics of real systems. 
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